Study of the thermal abelian monopoles with proper gauge fixing 
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The properties of the thermal abelian monopoles are studied in the deconfinement phase of the 
SU (2) gluodynamics. To remove effects of Gribov copies the simulated annealing algorithm is applied 
to fix the maximally abelian gauge. Computing the density of the thermal abelian monopoles in the 
temperature range between 1.5 T c and 6.9 T c we show, by comparison with earlier results, that the 
Gribov copies effects might be as high as 20% making proper gauge fixing mandatory. We find that 
in the infinite temperature limit the monopole density converges to its value in 3-dimensional theory. 
To study the interaction between monopoles we calculate the monopole-monopole and monopole- 
antimonopole correlators at different temperatures in the region (1.5 T c ,6.9 T c ). Using the result of 
this study we determine the screening mass, monopole-monopole coupling constant, monopole size 
and monopole mass. In addition we check the continuum limit of our results. 
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I. INTRODUCTION 

Our study of the thermal monopoles is motivated by 
the hypotheses that some of the quark-gluon plasma 
properties may be dominated by a magnetic compo- 
nent PHI], at least for the range of the temperature 
values close to the transition. The monopoles or mag- 
netic vortices might be responsible for the very low 
viscosity indicating that the quark-gluon plasma is an 
ideal liquid. 

The way one can study the monopoles properties 
in lattice nonabelian gauge theories is via an abelian 
projection after fixing the maximally abelian gauge 
(MAG) 0,1. Since first lattice studies of this gauge 
0, 0] the properties of the abelian monopole clus- 
ters had been investigated in a numerous papers both 
at zero and nonzero temperature. The evidence was 
found that the nonperturbative properties of the gluo- 
dynamics such as confinement, deconfining transition, 
chiral symmetry breaking, etc. are closely related to 
the abelian monopoles defined in MAG |8l— Tl0j| . This 
was called a monopole dominance. 

It was shown in Ref. @] that thermal monopoles 
in Minkowski space are associated with Euclidean 
monopole trajectories wrapped around the temper- 
ature direction of the Euclidean volume. So the 
density of the monopoles in the Minkowski space is 
given by the average of the absolute value of the 
monopole wrapping number. First numerical investi- 
gations of such wrapping trajectories were performed 
in Refs. [HI and [12] ■ A more systematic study of 
the thermal monopoles in 517(2) Yang-Mills th eory 
at high temperature has been performed in Ref. [13 1. 
It was found in [l3| that the density of thermal 
monopoles p was independent of the lattice spacing 



as it should be for a physical observable. It was con- 
cluded that the monopole density was well described 
by p oc T 3 /(logT/A e// ) 2 with A eff ~ 100 MeV while 
the behaviour p oc T 3 /(logT/A e //) 3 , predicted by di- 
mensional reduction arguments, was compatible with 
data for T > 5 T c . 

The density-density spatial correlation functions 
has been also computed in [l3j]. It has been shown 
that there is a repulsive (attractive) interaction for a 
monopole-monopole (monopole-antimonopole) pair, 
which at large distances might be described by a 
screened Coulomb potential with a screening length 
of the order of 0.1 fm. 

It is known that the Gribov copies effects are strong 
in the MAG [3]. In [l4| a conclusion has been 
made that results for gauge noninvariant observables 
can be substantially corrupted by inadequate gauge 
fixing. For nonzero temperature the effects of Gribov 
copies were not investigated so far. In this paper we 
would like to close this gap. Following [lj] we ap- 
ply the simulated annealing algorithm with 10 gauge 
copies for every configuration to solve the problem of 
Gribov copies by approaching the global maximum of 
the gauge fixing functional. We present our results 
for the density of the thermal monopoles and the cor- 
relation functions and compare them with the results 
of Ref. [13| . We show that the Gribov copies effects 
are indeed large. We also present results for screening 
mass and coupling constant of monopoles interaction, 
the size and mass of monopole. 

The paper is organized as follows. In the next 
section we briefly review the details of the simula- 
tion. In the section III we study the dependence of 
monopole density on the temperature. In section IV 
the monopole-monopole and monopole-antimonopole 
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correlators are calculated at different temperatures. 
Using the result of this calculation we determine the 
size of monopole, Debye screening mass, monopole- 
monopole coupling constant. In section V we calculate 
the dependence of the monopole mass on the temper- 
ature. In section VI we consider the continuum limit 
of our results. In last section the results of this paper 
are summarized. 



II. SIMULATION DETAILS 

We study the SU(2) lattice gauge theory with the 
standard Wilson action 



s = p E E 1 - \ ^{u x ^u x+ ^ui + ^ui 

x [i~>ii L 



where (3 = 4/<?g and go is a bare coupling constant. 
The link variables U xli G SU(2) transform under 
gauge transformations g x as follows: 

Uxv ^ U° M = giU x »g x+ ^ ■ g x G SU(2) . (1) 

Our calculations were performed on the asymmetric 
lattices with lattice volume V = L t L^, where L t s is 
the number of sites in the time (space) direction. The 
temperature T is given by 



T 



1 

aL t 



(2) 



where a is the lattice spacing. 

The MAG gauge is fixed by finding an extremum of 
the gauge functional 



Fu(g) 



(3) 



with respect to gauge transformations g x . We apply 
the simulated annealing (SA) algorithm which proved 
to be very efficient for this gauge [l4| as well as for 
other gauges such as center gauges and Landau gauge. 
To further decrease the Gribov copy effects we gener- 
ated 10 Gribov copies starting every time gauge fixing 
procedure from a randomly selected gauge copy of the 
original Monte Carlo configuration. 

In Table U we provide the information about the 
gauge field ensembles used in our study. 



III. THERMAL MONOPOLE DENSITY 

We locate the wrapped monopole trajectories us- 
ing a standard prescription. The monopole current 
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TABLE I: Values of /3, lattice sizes, temperatures, num- 
ber of measurements and number of gauge copies used 
throughout this paper. To fix the scale we take yfd = 440 
MeV. 



is defined on the links {a;, //}* of the dual lattice 
and take integer values j^ix) — 0, ±1,±2. The 
monopole currents form closed loops combined into 
clusters. Wrapped clusters are closed through the lat- 
tice boundary. For a given cluster the wrapping num- 



ber N wr is defined as 



A^ mr — 



- E 



3i(x) 



(4) 



j4(x)£cluster 



It takes values 0,±1,±2.... 

The density p of the thermal monopoles is defined 

as 



E 



clusters 



L?a 3 



(5) 



In Figure Q] we show our data for the set of lattices 
listed in the Table HI The data of Ref. [l3j are also 
shown for comparison. One can see strong Gribov 
copies effects from the Figure [TJ They are about 
20% at T/T c = 2 and increase to almost 30% for 
T/T c = 6.9. The dimensional reduction suggests for 
the density p the following temperature dependence 
at high enough temperature 



p(T) = {c P g\T)Tf 



(6) 



where the temperature dependent running coupling 
g 2 (T) is described at high temperature by the two- 
loop expression with the scale parameter A^: 



^ln(T/A T ) + i ^(ln(21n(T/A T )) (7) 
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In Ref. [13( the data for density p were fitted to a fit 
function 



(8) 



(log(T/A e// ))< 



3 




which is motivated by a one loop expression for g 2 (T) . 
A good fit was obtained for T/T c > 2 with a = 2 and 
it was noted that a = 3, which corresponds to the one 
loop expression for g 2 {T) is compatible with data for 
T > 5 T c . 

We fit our data to function (HJ),® with fitting pa- 
rameters c p and At- The reasonably good fit with 
X 2 /dof = 1.4 was obtained for T > 3T C with val- 
ues for fit parameters c p = 0.79(1), T C /A T = 2.52(6). 
Thus behaviour of the density p(T) is well described 
by the fitting function motivated by dimensional re- 
duction prediction. However, for the spatial string 
tension (T) the parameter At was found very dif- 
ferent T C /A T = 13.16(17). This implies that 
the dimensionless ratio p 1 / 3 {T) / a\^ 2 (^T) is decreasing 
with temperature. 1 This decreasing can be indeed 
seen from Figure [2] where this ratio is depicted. For 
al^ 2 (T) we took the fitting function obtained in 15]: 



^jT)=c s g 2 (T)T 



(9) 



with c s = 0.369(14) and T C /A T = 13.16(17). 

Note that decreasing seen in Figure [5] is slow, the 
ratio changes only about 10% over our range of tem- 
peratures. We fitted the ratio p 1/3 {T)/al /2 (T) to a 
polynomial fit 
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R + R 2 /x 2 , x = T/T c 



(10) 



for T/T c > 4 and obtained R = 0.49(1). This value 
to be compared with respective value for 3d SU{2). 



The monopole density p$ in this theory was computed 
in [13 : P3 = 0.0078(2) ff f. Let us note that this 
value was obtained with the overrelaxation gauge fix- 
ing algorithm which, as is well known from 4D studies 
14j . gives overestimated value for the density. Using 
the value for 3c? string tension er 3 from [Tj| y/as = 

0.3353g| we obtain for the ratio p\ /a\ = 0.59(1). 
This is in good agreement with our value for R if we 
take into account that the value for the 3<i density is 
overestimated as was discussed above. 

Thus in this section we have demonstrated that 
the density of monopoles changes with temperature 
in good agreement with predictions of the dimen- 

3/2 

sional reduction and being measured in units of ay 
is not far from its infinite temperature limit even at 
T/T c = 1.5. 



IV. STUDY OF THE INTERACTION 
BETWEEN MONOPOLES. 

A. The details of the calculation. 

In this section we study the interaction between 
wrapped monopoles, through the measurement of the 
spatial correlation function g(r) — (p(O)p(r)) / p 2 for 
the thermal monopole- (anti)monopole pairs. On the 
lattice this correlation function can be measured as 
the ratio of the number of the (anti)monopole located 
at distances (r, r + dr) from a given monopole and the 
number (anti)monopoles in the volume 4irr 2 dr in a 
completely homogeneous system 



7(r) 



1 dN(r) 
p Ai:r 2 dr ' 



(11) 



1 (Note that for magnetic screening mass T c /At = 5.1(9) was 
found in [3). 



where dN(r) is the number of particles in a spherical 
shell (r, r + dr), p is the average density of monopoles. 
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FIG. 3: The correlation functions g(r) for T/T c = 1.5 
(squares) and 4.8 (circles) for the monopole-monopole 
(empty symbols) and monopole-antimonopole (full sym- 
bols) cases. 

In order to take into account discretization effects we 
calculate the volume 4irr 2 dr as a number of three 
dimensional cubes on the lattice which are located 
in the spherical shell (r,r + dr). Note that at large 
distances the densities p(0),p(r) become uncorrelated 
and g(r) r ^oo -> 1 

We have measured the correlation functions at the 
following temperatures T/T c = 1.5,2.0,3.0,4.8,6.9 
where our statistics was large enough. The plots of the 
correlation functions at temperatures T/T c = 1.5,4.8 
are shown in Figure [3] It should be noted that the 
plots of the correlation functions shown in Figure [3] 
are in agreement with the results of paper [l3| where 
the same correlation functions were studied. 

The correlation functions g(r) contain information 
about interaction properties of monopoles in QGP. In 
order to get an idea about monopole-monopole and 
monopole-antimonopole interaction potentials we use 
the ansatz based on Boltzmann distribution 

ff(r)=e*p(-^), (12) 

where U(r) is the interaction potential of monopoles 
in QGP. Having the correlation function one can ex- 
tract the interaction potential U(r). The interaction 
potentials of monopole-(anti)monopole pairs for the 
lowest temperature T/T c = 1.5 and for the tempera- 
ture T/T c = 4.8 are shown in Figure|U The potentials 
of the monopole-(anti)monopole pairs interaction at 
other temperatures studied in this paper are similar 
to that at temperatures T/T c — 1.5,4.8. 

B. The size of monopole. 

From Figure |4] it is seen that monopole-monopole 
repel from each other for all values of r. The 
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FIG. 4: The potentials U{r) for T/T c = 1.5 and 4.8. The 
symbols are the same as in Fig. [3] 



interaction potential of monopole-antimonopole has 
rather complicated structure. At large distance the 
monopole-antimonopole attracts. However, the po- 
tential reaches minimum and after the minimum it 
rises, becoming repulsive. 

Qualitative features of the monopole-antimonopole 
potential are similar to that in dilute one atomic gas. 
In the later case the interaction between molecules 
is attractive at large distances, reaches minimum at 
some point and becomes repulsive at distances of or- 
der of the size of molecule. From this analogue one 
can assume that monopoles have some size and inter- 
actions becomes repulsive at distances of order of the 
monopole size r mon . In this paper we determine the 
size of the monopole as 2r mon = r m i n , where r m i n is 
the position of the minimum of the interaction poten- 
tial. In Figure [5] we plot the product r mon ^fal as a 
function of temperature. Is is seen that within the er- 
ror of the calculation the plot does not contradict to 
the behaviour which is suggested by the dimensional 
reduction at large temperatures. However, because of 
the uncertainty we cannot state that the dimensional 
reduction really takes place. To summarize the result 
of this study: monopole is not a pointlike object but 
an object with nontrivial core with the size ~ 0.05—0.1 
fm depending on the temperature. 



C. The tails of the correlation functions. 

Now let us pay attention to the tail of the correla- 
tion functions at r — > oo. In the numerical analysis 
we are going to fit the tail of the correlation function 
as it is fitted in the electromagnetic plasma 

g(r)=exp(-^e- m ^, (13) 
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FIG. 5: The product r m0 n^fol as a function of tempera- 
ture. 



where vhd is the Debye screening mass, a is the 
monopole coupling constant. Now two comments are 
in order 

• In electromagnetic plasma the constant a is the 
square of ion charge. However, in QGP the in- 
teraction is much more involved. In particular, 
the mechanism of Debye screening is not the 
same as it is in electromagnetic plasma. What 
is more important: in view of our recent find- 
ing about dyon nature of monopoles [l9|, one 
can state that monopoles interact both by chro- 
moelectric and chromomagnetic charges. These 
facts mean that the constant a cannot be con- 
sidered as a square of chromomagnetic charge. 

• As was noted above correlation functions (|TT1) 
are normalized as 5(r) r ^. 00 — > 1. However, 
in lattice simulation the correlation functions 
at large distances differ from unity by amount 
~ 0.01. Note the same effect was observed in 
paper [20]. Analysing available data we have 
come to the conclusion that this is finite volume 
effect. In order to take into account this effect 
we fit our data by the function 



g(r) = exp 



Tr 



(14) 



It is natural to expect that in the limit r — > oo 
monopole-monopole and monopole-antimonopole cor- 
relators have equal Md and up to the sign equal 
a. Unfortunately we cannot check this since our lat- 
tices have finite sizes. It is also important that the 
monopole-antimonopole correlator has rather nontriv- 
ial behaviour changing from attraction to repultion at 
some point r m i n . Evidently to get reliable estima- 
tion of the Md and a at the tail of the correlator one 
must consider the distances r 3> r m i n what is also 



impossible due to the finite lattice sizes. For this rea- 
son we determine the parameters Md and a from the 
monopole-monopole correlator. 

In Table |TT] we present the result of the fit for the 
monopole-monopole potential. It is seen that the fit 
is rather good. In Figure [5] we plot the Md/ ' ^fa~ s as a 
function oiT/T c . From this plot one sees that within 
the error of the calculation the ratio Mul yfol does 
not depend on the temperature what is suggested by 
the dimensional reduction at large temperatures. 

The authors of paper [20J studied the temperature 
dependence of the constant a. They found that the 
constant a rises with temperature. Our results con- 
firm the fact of rising the a. It should also be noted 
that the values of the constant a obtained in [2(| and 
in this work are a little bit different but in a reasonable 
agreement with each other. 

The authors of paper [2l[ found the value 
M-r>l\fo~ s = 1-77. As it is seen from Figure H] our 
result for this ratio is MpL/a^ ~ 2 what is in good 
agreement with paper [2l| . taking into the account 
the fact that actually the definitions of the Md in the 
both papers are different. In paper [l3| the Debye 
mass was estimated 1/Md ~ 0.1 fm. Our result for 
the 1/Md ~ 0.1 — 0.2 fm depending on the temper- 
ature, what is also in reasonable agreement with the 
estimation of paper (T3 |. 
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TABLE II: The masses Md and coupling constants a ob- 
tained by fitting of the correlators to eq. (|13p . 



V. THE MONOPOLE MASS. 

There are a lot of definitions of the monopole mass 
(see [22]). In this paper we are going to apply the def- 
inition which is based on the formalism of functional 
integration. Consider the trajectory of free nonrela- 
tivistic particle of mass m. The spatial fluctuation of 
the trajectory of such particle 



Ar 2 = T 



f 1/T dt((r(t)-r(0)) 2 ), (15) 
Jo 



is directly connected to the particle mass 
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FIG. 6: The ratio of Md to ^fal as a function of temper- 
ature. 

The lattice version of the equation (fT5|) can be written 

as 

Ar 2 = ^« f .^i) 2 )- ( 17 ) 

Where the sum is taken along the monopole trajec- 
tory, L is the total length of the trajectory. 

With equations ([To]) . (JTTj) we get the result pre- 
sented in Table IIIII Kinetic energy of monopole is 
~ T. Monopole mass is ~ 2. IT at T/T c = 1.5 reach- 
ing ~ 8.2T at T/T c = 6.9 what legitimates nonrela- 
tivistic approximation. It should be noted that the 
authors of paper [22| also calculated the mass (p~6|) . 
Their result approximately 30% larger than ours. 
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TABLE III: The monopole masses m and spatial fluctua- 
tions at different temperatures. 

It should be noted that the larger the temperature 
the more static monopole trajectories. What means 
that contrary to the behaviour of the screening mass 
Mb/T the monopole mass m/T increases with tem- 
perature. In Fig. [7] we plot the monopole mass as 
a function of temperature T/T c in logarithmic scale. 
From this figure one sees that all points beautifully 2 
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FIG. 7: The monopole mass m/T as a function of temper- 
ature T/Tc in a logarithmic scale. 



lie on the line 

with b = 3.653(6), A m /T c = 0.718(2), x 2 /dof ~ 0.2. 
This means that m/T ~ l/g 2 {T) or Ax 2 T 2 ~ g 2 (T). 



VI. THE CONTINUUM LIMIT. 

To check the continuum limit of our results we have 
generated 1000 configurations of lattice 8x48 3 with 
(3 = 2.635. The temperature of these configurations is 
T/T c — 1.5 and it is equal to the temperature of the 
configurations at lattice 4x36 3 with f3 = 2 A3 but the 
lattice spacing is two times smaller. 

First we note that the monopole density calculated 
at lattices with a = 0.1 fm and a — 0.05 fm differs 
from each other by ~ 5%. What can be considered 
as a good continuum scaling and that even at lattice 
with a = 0.1 fm we are not far from the continuum 
limit. 

Further let us consider monopole correlators. In 
Fig. |S]we plot the correlators at different lattice spac- 
ings. From this plot one sees that both monopole- 
monopole and monopole-antimonopole correlators are 
in good agreement with each other at distances rT > 
0.7. At distances rT < 0.7 the correlators differ 
from each other. However, this difference can be at- 
tributed to the fact that the lattice spacing a = 0.1 
fm is too large and the lattice to rude to probe in- 
ternal monopole structure. At the same time the lat- 
tice spacing a — 0.05 fm is small enough to do this. 



2 Note that the uncertenties of the data shown in Table HITl are 



very small. 
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FIG. 8: Monopole-monopole and monopole-antimonopole 
correlators at temperature T/T c = 1.5 for lattices 8x48 3 
with a = 0.1 fm and /3 = 2.635 with a = 0.05 fm. 



copies the simulated annealing algorithm is applied 
for fixing the maximally abelian gauge. 

We calculated the density of thermal abelian 
monopoles at different temperatures. It was found 
that at temperatures T > 3T C our data can be very 
well fitted by the function which is motivated by di- 
mensional reduction prediction p ~ (g 2 (T)T) 3 . 

I addition to the monopole density in this paper we 
studied the monopole-(anti)monopole spatial correla- 
tion functions. We calculated the correlation func- 
tions at the temperatures T/T c = 1.5, 2.0, 3.0, 4.8, 6.9. 
This result allowed us to determine the potentials of 
monopole-monopole and monopole-antimonopole in- 
teractions. As one can expect monopoles repel from 
each other. The tails of the interaction potential can 
be well fitted by the screened Coulumb potential 



(20) 



Note that the domain of repultion of the monopole- 
antimonopole pair is clearly seen on this plot. 

In section IV the size of monopole was defined as 
2r mo „ = ro, where rg is the position of the maxi- 
mum in the monopole-antimonopole correlator. The 
position of the maximum shifted to the right but not 
greater than by 20%. If the maximum were lattice 
artefact for dimensional arguments its position would 
shift to the left. From this we conclude that monopole 
has some physical size which has good scaling be- 
haviour. 

The tail parameters a,Afo/T ([M)) and monopole 
mass m/T ([To]) at lattice with a = 0.05 fm have the 
following values 



a = 3.0 ±0.8, 



Md 
T 



1.9 ±0.2, —=3.96 ±0.2. 
T 

(19) 

We see that taking into the account uncertainties of 
the calculation the results for the a^Mu/T are in 
agreement with Table flTJ As to the monopole mass 
m/T its value is by 30% larger than the result pre- 
sented in Table IIIII Note that Table IIIII was obtained 
at lattices with N t = 4, but results Table [JJJ] were ob- 
tained at lattice with N t — 8. Evidently, the thermal 
monopole trajectory parameters are very sensitive to 
the Nt. So, this could be the reason of the deviation 
of the both results. 



VII. CONCLUSION. 

In this paper the properties of the thermal abelian 
monopoles are studied in the deconfinemcnt phase of 
the SU(2) gluodynamics. To remove effects of Gribov 



We determined the dependence of the parameters 
a, Mjj on temperature. 

The interaction potential of monopole- 
antimonopole has rather complicated structure. 
At large distances it is attractive. However, the 
potential has a minimum at some distance and after 
the minimum it rises, becoming repulsive. The 
position of the minimum can be considered as a 
double monopole size. In this way we estimated the 
size of monopole which turned out to be 0.05 — 0.1 
fm depending on the temperature. 

We have also calculated the mass of monopole and 
the dependence of the mass on temperature. 

The last point considered in this paper is the con- 
tinuum limit of the results obtained in this paper. The 
calculation was done at lattice 8x48 3 with two times 
smaller lattice spacing a — 0.05 fm in comparison to 
the lattice 4x36 3 (see Table [TJ. All results except 
the monopole mass have good scaling behaviour. The 
monopole mass is increased by 30%. To undestand 
the origin of this increase further study is needed. 
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